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1 Introduction 



In the recent days, our perspectives in the search for physics beyond the Standard Model (SM) have enriched 
considerably with the developments of theories with extra (compact) dimensions [1-3]. In particular, these extra 
spatial dimensions allow to lower the gravity scale down to the TeV [1,2], thus addressing the long standing 
puzzle of the gauge hierarchy problem, when gravitational interactions feel these dimensions. Such TeV-scale 
Extra Dimensional (ED) models have attracted a great attention from the experimentalist community, since 
one can undertake their precision studies in the next generation colliders. Depending on the geometry of the 
extra dimensions (shape and/or size), and the field localization, a wide variety of extra dimensional models have 
been proposed. In addition to the gauge hierarchy problem, these models address several other questions in the 
different avenues of particle physics, cosmology and string theory. 

Among a large number of possible realizations of the ED models, the scenario, with additional right handed 
neutrinos (being gauge singlet, and thus sterile) in flat compact extra dimensions and SM fields localized on a 
brane, provides an interesting alternative [4-7] capable of accounting for the observed light neutrino masses [8] . 
Also, the detailed experimental results on neutrino properties, in the context of this scenario, may be utilized 
to shed light on the bulk geometry (for phenomenological studies, see [4, 5,9]). 

In this paper we study this scenario based on the formalism of [11]. Our first aim is to determine whether a 

realistic neutrino spectrum can still be gcnierated when the constraints on mixing angles between active and 
sterile neutrinos are taken into account, and particularly the severe constraint from the SNO experiment. Such 
constraints should indeed apply on the considered scenario, since the new right handed neutrino, as well as each 
component of its infinite massive tower of Kaluza-Klein (KK) excitations, constitute sterile neutrinos which mix 
with SM active neutrinos. We consider the general Lagrangian, which possess two distinct sources of Lepton 
number (L) breaking. The phenomenological constraints on neutrinos were applied, in [12], to the model with 
an L-conserving Lagrangian assuming no brane shifting. 

Our second goal is to address the issue of a successful leptogenesis within such realistic ED model in agreement 
with experimental neutrino data. More precisely, by studying the decays of the heavy KK neutrinos, we 
explore the possibilities that these scenarios account for the observed baryonic asymmetry of the Universe by 
means of the Fukugita-Yanagida mechanism of leptogenesis [13]. In this mechanism an excess of L, generated 
by out-of-equilibrium i-violating decays of the heavy neutrinos (arising with the i-violating terms present 
in our Lagrangian), is converted into a Baryon rnimber (B) asymmetry through (B + L)-violating sphaleron 
interactions. At this stage we mention a preliminary work [6] where the leptogenesis has also been considered in 
the same context of flat extra dimensions, but with an extended Higgs sector compared to the SM and without 
considering the various experimental neutrino data. 

We will consider the case where gravity may propagate in an higher dimensional space of [l + (3+ng)] dimensions, 
where rig > 1. Hence we have two variants of the scenario, one with Ug = 1 and the other with ng > 1. In all 
cases, the SM particles are localized on a [1 + 3] dimensional subspacc (3-brane). For simplification reason, the 
sterile neutrinos are assumed to propagate in [1 + (3 + 1)] dimensions. 

The organization of the paper is as follows. We begin with a brief outline of the essential features of the extra 
dimensional scenario in Section 2. The relevant phenomenological constraints on neutrinos are discussed in 
Section 3, while our numerical results are presented and discussed in Section 4. Next in Section 5, the issue of 
leptogenesis is addressed. Lastly, we summarize our studies, and conclude, in Section 6. 
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2 Minimal Higher Dimensional Model for Neutrinos 



In this section we describe the minimal ED framework with a 5-dimensional iso-singlct right handed neutrino 
added to the field content of the SM [4, 11]. Here, we stick to one generation of SM neutrinos. The fifth flat 
dimension, along which propagates the right handed neutrino, is compactified over a jZi orbifold. The SM 
fields are localized on a 3-brane, whereas gravity propagates in the bulk. 
The leptonic field content is. 



L{x) 



Ir{x) , 



as in the SM, where v^, li,, are 4-dimensional Weyl spinors, plus the 5-dimensional singlet neutrino. 



N{x, y) 



fj{x,y) 



(1) 



(2) 



where y parameterizes the compact fifth direction, and ^, rj are the 2-component spinors. The y-coordinatc is 
compactified on a circle of radius R (the periodic boundary condition N{x, y) = N{x, y + 2ttR) is imposed) on 
the singlet neutrino field. Additionally, under the action of Z2 symmetry, the two 2-component spinors may be 
associated with opposite parities, 



£.{x,y) = ^{x,-y), vix,y) = -T}{x,-y). 



(3) 



The SM fields are localized at the brane, which, just for generality, we may assume to be at y = a, instead of 
at the orbifold fixed point y = 0. For this model, an effective Lagrangian respecting the 4-dimensional Lorentz 
invar iance reads as, 



2-kR 



jCeS = jdy |iv(n''5^ + -i^dy^N - i (mN^C^^'^-^-i^N + h.c.) 


^1 ri*. , h2 



+ 5{y-a) 



(Mir)V2 



L^*r] + h.c. 



6{y- a)£sM 



(4) 



where $ = ia2^* is the hypercharge-conjugate of the SM Higgs doublet with hypercharge = 1, £sm 

is the SM Lagrangian, M is the Majorana mass for A'' (we do not specify its scale for the moment), C^^-' is the 
5-dimensional charge conjugation operator and Mp is the fundamental rig-dimensional gravity scale, given by 

Mp = {2-kMpRYoI^Mf, (5) 

for the simple case where all the compactification radii are of equal size R, Mp being the 4-dimensional Planck 
scale. A Dirac mass term moNN is not allowed in Eq.(4) because of the Z2 discrete symmetry. Note that the 
Majorana mass term for A^ breaks the 5-dimensional Lorentz invariance as in [10]. 
Following Eq.(3), the 2-component spinors ^ and r] can be expanded in Fourier series as. 



^{x,y) = 



I-kR 



V2TrR 

V{x,y) = ^ 77„(a;) sin 



cos 



ny 
R 



(6) 



(7) 



where the chiral spinors ^n{x) and r]n{x) form an infinite tower of KK fields. Using these expansions and 
integrating out the ^/-coordinate, the effective Lagrangian reduces to. 



E 



(8) 
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in a basis in which M is positive, and with: 
{n) _ hi 



h 



(n) 



(27rMFi?)i/2 



sm 



na ^ 




( Mf 


li) 




V Mp 


na \ 




( Mf 


li) 







l/ng 



l/Ug 



hi cos 



/i2 sin 



na 
na 



= hi cos 
= /i2sin 



na 
na 



(9) 
(10) 



For deriving the last two equahties on the right hand sides of Eqs.(9)-(10), we have made use of Eq.(5). 



Eqs.(9) and Eq.(lO) tell us that the reduced 4-dimensional Yukawa coTiplings h'i^l can be suppressed by many 
orders depending on the hierarchy between Mp and Mf; for example, if gravity and the bulk neutrino feel the 
same number of extra dimensions, say Ug = 1, then these couplings are suppressed by a factor Mf/Mp 10~^^, 
for Mf « 10 TeV (see also [4, 5]). 

From Eq.(2) we observe that ^ and fj have the same lepton number. Thus, the simultaneous presence of the two 
operators L$*^ and L^*r] in Eq.(8) leads to lepton number violation. If the brane were located at one of the 
two orbifold fixed points (y = or wR), the operator L^*r] would have been absent in Eq.(8) as a consequence 
of the discrete Z2 symmetry. The two operators can coexist, leading to breaking of the lepton number, if we 
allow the brane to be shifted by an amount a{y^ 0) from the orbifold fixed points. In fact, it is possible to 
perform such a shifting of the brane (even in a continuous way) respecting the Z2 invariance of the original 
higher dimensional Lagrangian under certain restrictions in Type-I string theories [14]. As indicated in [4, 11], 
the Z2 invariance can be taken care of by allowing the replacements. 



^Siy-a) 
r]S{y-a) 



i C 6{y-a) + 6{y + a- 2t,R) 



2^ 



5{y -a) - 6{y + a- 2-kR) 



(11) 



with < a < 7ri? and < 2/ < 2ttR. Obviously, a Z2-invariant implementation of brane-shifted couplings would 
require the existence of at least two branes placed at y = a and y = 211 R — a. A remarkable feature of the 
brane-shifted framework was pointed out [11], where it has been shown that in such a framework it is possible 
to completely de-correlate the effective Majorana- neutrino mass (to) , that is effectively measured, and the scale 
of light neutrino masses, as to have (to) within an observable range. 

Therefore, the Lagrangian (8) contains two types of Majorana neutrino mass term (involving respectively the 
parameters M and h^'') which lead both to a breaking of L. This L-breaking is a necessary ingredient of 
leptogenesis. 

Following the notations of reference [4], we now introduce the weak basis for the KK Weyl-spinors, 

1 



X±n = (Cn ± rtn), 



(12) 



followed by a rearrangement of ^0 and Xn states, such that, for a given value of k (say, k = fco) the smallest 
diagonal entry of the neutrino mass matrix is 



£| =min( M - ^ ) < l/(2i?) 



In this new basis the reordered 4-component Majorana-spinor vector is the following. 



Xvi 



Xko 
Xko 



Xko + 1 
Xfeo-l-1 



Xfeo-1 
Xfeo-1 



Xko+n 
Xko+n 



Xko- 
Xko- 



while the effective Lagrangian for right handed neutrinos reduces to. 



(13) 



(14) 



(15) 
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where A4^^ is the corresponding neutrino mass matrix given by, 



( ^ 






1) 






2) 


m(2) 




e 


















m(-i) 





£ — 
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R 













m(i) 








£ + i 










m(-2) 













£ — 


2 
R 





m(2) 


















^+1 



V 



(16) 



The most important consequence of such a rearrangement is that the mass scale M, which we did not specify 
earUer but which could be as large as possible, is now replaced by the light mass scale e. The mass entries in 
the first row and first column of matrix (16) are given by the relation. 



m 



(») - 



V 



hi cos 



(n — fco)a 
R 



sin 



{n — ko)a 
R 



with, 



m cos 



na 
'R 



(17) 



m 



= tan 




where v is the vacuum expectation value of the SM Higgs boson. 



(18) 
(19) 



We observe that the structure of neutrino mass matrix (16), originating from a combination of Major ana and 
Dirac mass terms, is different from the pure Dirac mass case considered in the preliminary work [12]. 



3 Constraints for a Realistic Spectrum 

The neutrino spectra is obtained by diagonalizing the mass matrix A4^^ of Eq.(16). It has been shown in [11], 
for a one-generation case, that the eigenvalue equation for a restricted class of cases with a = nR/q, where q is 
an integer greater than 1, is easily tractable analytically. Wc have, in the generic case, solved the characteristic 
eigenvalue equation numerically to obtain the eigenvalues of M^^, sticking to the one-generation case. We have 
diagonalized the N x N matrix of Eq.(16) for a N satisiying Mp > Mmax (see end of the paragraph). We have 
systematically chocked that the N value used (typically ^ 500 for the smallest R~^ values to ^ 10 for the largest 
R~^ values) constitutes a dimension at which the lightest eigenvalues converge and are completely stabilized 
(with increasing TV). We have performed a scan over the whole parameter space of the ED model described 
above in order to find the regions in agreement with basic experimental c;onstraints on neutrinos. The total 
parameter space consists of R^^, £, the complex Yukawa couplings hi and /i2, the brane-shift parameter a, and 
the phase (ph defined in Eq.(19). From Eq.(5), the effective cut off scale Mp for a given R~^ is completely fixed 
for Ug = 1; while for rig > 1, it depends on Ug and the respc;c;tive compactification radii (respecting Mp > Mmaxj 
where Mmax is the mass of the heaviest neutrino eigenstate considered). 

While scanning over the total parameter space described above in order to search for a realistic neutrino mass 
spectrum, we have to ensure several constraints described in the following; 

• There exist phenomenological constraints on the mixing angles between active neutrinos and sterile ones. 
These constraints apply on the present model, as the 0-mode and KK excitations of the additional right 
handed neutrino behave exactly as sterile neutrinos — the right handed neutrino having no electroweak 
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interactions (see e.g. [15]). 

For instance, in a model with a unique sterile neutrino, and assuming for simplicity one lepton fla- 
vor, the mixing angle 9s between the active and sterile neutrino is constrained typically by tan^^ < 
10"'' - 10-1 for Am^ e [lO-^^ - 10^] eV^ (Am^ being the mass eigenvalue squared difference) [16] from 
cosmological and astrophysical considerations combined with data from atmospheric, solar (including 
SNO), reactor and short base-line experiments. In the limiting case Am^ ^ Am^sun ^ 10"^ eV^, the SNO 
bound on the fraction of oscillating sterile neutrinos coming from the sun = siii^ 9s/2 [16], becomes 
approximatively [17, 18] 

ris < 1.2 X 10-1 at la. (20) 

Similarly, in our framework, one can also take as an approximation that the Xka state, whose mass e is the 
first appearing as a diagonal matrix element in Eq.(16), is the only sterile neutrino that may not decouple 
from the active SM neutrino f^. Indeed, by definition of e in Eq.(13), the following Xko±n states get larger 
and larger diagonal mass matrix elements e ±n/R and thus decouple more and more from vg, given the 
symmetric structure of the mass matrix in Eq.(16). Therefore, under the above approximation, the SNO 
bound on Og apply in our framework also, where the mixing angle 9s is deduced only from the 2x2 block 
matrix of Eq.(16): 

tan26»sRi— (21) 

The above SNO bounds impose a significant upper limit on 9s leading to a case comparable to the see-saw 
mechanism [19] (see also the third point) so that the smallest eigenvalue m^i and the next-to-smallest 
eigenvalue rriiy^ read as 

TOi/j « m'-'^', TOi^2 ~ £. (22) 

The constraint in Eq.(20) imposes m^^^ <C e {c.f. Eq.(21)), which may be achieved in the following two 
cases (see Eq.(17)): 

m<C£, and/or i^/, ~ 7r/2 -)- 2g7r, (23) 

where q is an integer. 

The previous approximations regarding the implications of the constraints on sterile neutrino mixing will 
turn out to be useful in the interpretation of our results. In order to quantitatively take into account 
the effects of the KK tower of neutrino states when applying the main SNO constraint, we have imposed 
the bound: P{i'e '-'s) < 0.40 at 90% C.L. (coming from the SNO data [20] combined with the Super- 
Kamiokande results [21]). P(fe — > I's) represents the transition probability between solar SM electron- 
neutrino Ve and sterile neutrino i^s, and, is given by P(z/(, ;/,) = 1 — J2e=e /i t -^('^<^ ^ where 
P{ve t^e) = J2n l^enl^jV^nl^ [22] with ^ = 6 in the one active flavor case. The sum here is taken over the 
index n, which labels the neutrino mass eigenstates, and, V diagonalizes the neutrino mass matrix Eq. 
(16). The sterile KK excitations are little affected by matter effects [12]. Following the same numerical 
approach for N as described in the beginning of this section, we have checked that for the chosen value of 
N (corresponding to the maximum value of n truncating the sum involved in the P(fe ^i) expression), 
P{ve ^ i^s) reaches a value which is not significantly sensitive to further increases of the mass matrix 
dimension. In other terms, a satisfactory level of convergence can be reached in the summation of KK 
excitations. In fact, the stabilized values of P{ue I's) for our case turn out to be systematically several 
orders of magnitude smaller than the combined SNO-Supcr-Kamiokande bound mentioned previously, in 
the final regions of parameter space which respect all the constraints discussed in this section. 

Due to the see-saw structure of the mass matrix, the two eigenstates ui and 1/2 are respectively composed 
mainly by the SM neutrino and sterile state Xfeo • Hence one can apply the phenomenological bound on 
nivi derived for SM neutrinos, which in the case of one flavor is typically [23] 



10(-2 - 0) gy_ (24) 



6 



• We have also ensured that the mass of the next-to-hghtest^ eigenstate V2 is larger than the mass of the 
Higgs boson, typically 

> 200 GeV, (25) 

for our purpose which deals with its possible Icpton-number violating decay into Higgs and leptons: 
U2 £^(jy^ (the three-body decay channel of the is 1^2 £"^£^1^1). 

The choice of this condition m^^ > 200 GeV is important for our study of leptogenesis (in Section 5). We 
could also have imposed this constraint on any higher-mass neutrino eigenstate jUi,. > 200 GeV with i > 2, 
and studied the decay i^i l^cjy^ . However, any higher state t-j (i > 2), in addition to its decay channel 
Vi — > possesses several other decay channels: their decays into lighter eigenstates Vi — > Vjl'^£~ or 

Vi VjViVi with i > j > 2. How the presence of these additional decay channels is going to affect the 
Boltzmann equations and the washout factors, and thus the resulting lepton abundance, is a matter of 
non-trivial calculation that we will not address here. Such additional decay channels do not exist for the 
state 1^2, as the decays 1^2 ^ i^2i^ii^i and 1^2 ^ 1^2^^^" arc not possible kinematically. Here we consider 
the possibility of the ED model to yield a successful leptogenesis via the 1^2 decay. 

With this bound m^^ ^ 200 GeV we are in the case where Am? = — mj^ » Aml^,^ so that the 
relevant SNO bound is the one given in Eq.(20). Furthermore, this condition forces mi,^ to be much larger 
than rriui and so constitutes another justification for approximating the two lightest eigenvalues by the 
see-saw formulas (22). 

At this stage the following remarks on neutrino mass are in order. There are typically three different mass 
suppression mechanisms in the extra dimensional scenario, that may suppress the mass of the SM neutrino. 

1. The usual see-saw like suppression, from a heavy Majorana mass M for the right handed neutrino. 

2. An effective ED see-saw type mechanism resulting from the massive KK tower [4,6], where the KK 
excitations of the additional right handed neutrino play the role of the right handed neutrino in the usual 

see-saw mechanism. 

3. A suppression of the Yukawa couplings from the ED wavefunction overlap mechanism affecting the neutrino 
mass via the ED volume factor {MpR}^^^ (see Eq.(18)) and the factor cos{na/R — <f)h) (see Eq.(17)). 

The main motivation for working in an ED scenario is to generate small SM neutrino masses using the second 
and third mechanisms only, i.e. the usual see-saw mechanism does not constitute the main suppression origin. 
Is that possible? As we have already discussed, the diagonal matrix element e contributes dominantly to the 
eigenvalue m,,^-^ since it is the lightest one. Hence in order to insure that the neutrino mass suppression is 
not mainly due to the usual see-saw mechanism, one can impose the condition M ^ e, which is equivalent to 
M > 1/2R (see upper left 2x2 block matrix in Eq.(16)). This condition can always be satisfied since the mass 
scale M never appears in this scenario directly: the presence of the infinite tower of KK states guarantees that 
only the value of M modulo i?^^, namely e (in absolute values) plays a physical role [4, 11]. 

We now discuss the ranges considered for each fundamental free parameter. First we made sure that it is mainly 
the second and/or the third mechanism (described previously) that suppresses m^^, by letting the numerical 
values of the Yukawa couplings hi and ft.2 to be within the natural range ~ 0.1 — 5 with the respective phases 
ranging between and tt. The phase (j)h, was also varied within and tt. We varied from a few tens of 
eV up to ~ 10^® GeV. For > 10^^ GeV, the heaviest state considered has a mass Mmax greater than Mp, 
i.e. our low energy eS'ective theory looses its predictivity. The range of R~^ being so huge, we considered a 
logarithmic binning of the scale of R^^. For each R~^, the parameter e may also vary within a large range, 
especially when R~^ is large, following Eq.(13). So we also scan over e by choosing randomly the values of a 
in £ = 10". The brane shift parameter a was scanned in a similar way in < a < ttR. Finally, we performed a 
scan by randomly choosing ~ 10^ points, so that each point corresponds to a distinct but random combination 
of all the parameters. 

^In our notations we order the mass eigenvalues increasingly, i.e. rrii/j < nii/j < ■ ■ ■. 
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4 Results and Discussion 



When scanning over the parameter space {R~^,e,(j)h,hi,h2,a} we have imposed a reahstic spectra, i.e. one 
respecting all the constraints mentioned in the previous section, by using the numerical estimation of this 
neutrino mass spectra. We have explored the two typical situations allowing for such realistic spectra (c./. 
Eq.(23)) within two scenarios characterized by ng = 1 and Ug > 1. Our observations and discussion for both 
the scenarios are given below. We have also included, in each case, possible analytical interpretations of our 
results. 



4.1 Scenario with rig = 1 

1. The scale Mp'- We recall (see Eq.(5)) that the cut off scale Mp depends only on for Ug = 1. As we 

increase the scale of from a few tens of eV to ^ 10^^ GeV, the corresponding scale M p increases from 
~ 10^ GeV to ~ 10^^ GeV, following Eq.(5). This makes the scale m in Eq.(18) increasing from ~ 10^ eV 
to ~ 10=^ GeV. 

2. The parameter R~^: We observed that the energy scale R~^ has to be > 400 GeV (corresponding to 

Mp > 10^^ GeV) for generating a realistic neutrino spectrum (see Fig.l). The reason is that R'~^ > 2 |e| 
(from Eq.(13)) with e ~ rrii.^ from Eq.(22) and we have imposed 712^2 > 200 GeV. With regard to the 
neutrino mass, what is the physical meaning of this lower bound on R~^7 We know from Eq.(5), that 
as R,^^ increases, Mp also increases but at a much slower rate compared to R~^. As a result, the 
quantity 1/y/MpR increases with R~^ so that the volume suppression mechanism for Yukawa couplings 
becomes less and less effective (see Eqs.(9), (10)), for a fixed coscph- Besides, the hierarchy between 
m(°) cx 1/^/MpR and R~^ increases with R^^, so that the effective ED see-saw mechanism becomes more 
effective in the reduction of m^^ if s is set at its maximum value 1 /2R. Indeed the denominator of m^^ 
(Eq.(22)) clearly increases at a faster rate compared to the numerator m^°^ - thus making the see-saw 
formulae more effective as R~^ increases. Therefore, the fact that there is a minimum value for R~^ means 
that the see-saw mechanism plays at least a certain role in the neutrino mass suppression. 

3. The pcirameter e: For all R~^ > 400 GeV, there exists a minimum value of e, starting from which 
we end up with a realistic spectrum for all possible higher values of e (< 1/2R). This minimum value is 
significantly larger than 200 GeV (recall that e « mi,^ > 200 GeV) for > 10^ GeV, as we see on Fig.l 
which shows the plot of s versus 1/R. This can be understood as follows. 

This figure results from the general scan, and thus includes both the cases m <S e and (j)h ^ tt/2 + 2q7r. 
The three straight lines correspond to the e values obtained analytically from the approximate see-saw 
relation in Eq.(22), for m^^ = 1 eV, 0.1 eV and 0.01 eV, and yTif+Ttf 4 and = (i.e. cos^/^ = 1). 
Therefore, the points obtained close to these lines correspond to the case where the m^^ suppression is 
mainly due to the configuration (1). The regions far below these lines indicate the simultaneous occurrence 
of both configuration (1) and (2) i.e. m -C e as well as ~ n/2 + 2q-K. Starting from these regions and 
approaching higher R~^ values, m increases, making the suppression from (1) less effective. Thus if 
we are to generate the same amount of suppression in m^j (or equivalently the same smallness of ratio 
nS^^ /e) for a fixed e, the TOyi suppression coming from case (2) must be more effective, i.e. has to be 
closer to 7r/2. This condition on 4>h basically reduce the allowed range for cjjh, so that the scan misses the 
allowed 4>h values. This explains the large vacant area in Fig.l below the three lines. Therefore, the easier 
way to generate the suppression in the higher R~^ region is by case (1), i.e. by pushing the minimum 
value of £ to be still higher to account for the increase in m. Quite naturally therefore, the available 
parameter region for higher R^^ decreases, and is concentrated near maximal values of e. Nevertheless, 
because of choosing to scan over the (j)h angle, the particular regions where qJcq -C R and h\ <C /i2 (leading 
to a 4>h value close to 7r/2) have been missed. Anyway such a hierarchy between the Yukawa couplings 
hi,h2 would have no theoretic;al explanation and, in c;ontrast, the; interest of the considered ED model is to 
have /ii ~ /i2 ~ 1 while suppressing the neutrino masses mainly from the second and/or third mechanism 
described above (the purely ED mechanisms). 
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Figure 1: Plot of e (in GeV) versus 1/R (in GeV) obtained from the scan. The three straight hues represent 
the analytical approximation of Eq.(22) for the indicated values of m^^. 

One must realize that the specific points in Fig.l, which are such that e <C correspond to a value 

of M extremely close to Uq/R (see Eq.(13)) and in turn to a situation of fine-tuning between the three 
parameters Af , fco a-nd R. 

4. The parameters (f)h and a: First let us discuss the spectrum. Since we are always working in a 
parameter region such that the off-diagonal elements of the mass matrix of Eq.(16) are much smaller than 
the diagonal elements, the eigenvalues m^^ of the matrix, up to a good approximation (see Sec A), are 
equal to the diagonal elements for i > 2: rriy^ ~ e, m^^ e — 1/ R, m^_^^ ^ e + 1/-R, and so on. In general, 
the masses of the higher neutrino states can be expressed as "rn^^^^^^^ ~ e — p/R, ra^^^^^^^ ~ £ + p/R, 
where p is an integer > 1. This clearly tells us that for a given R~^, the nature of the resulting spectra is 
primarily dictated by the choice of e. For example, if £ = l/2i?, the resulting spectrum becomes almost 
pairwise degenerate with \mi^^ \ ~ \fn^^\ (see Eq.(22)), |m^^| ~ Iwysl, and so on, while, for all the other 
choices of e, the resulting spectrum is hierarchical with m^^ < "^i^a < '^i'4 < '^1^5 < ■ • ■ 
In particular, we found that for a given R^^, the overall nature of the resulting spectrum feebly depends 
on the choice of the parameters (f>h and a. This is expected since these parameters appear only in the 
non-diagonal mass-matrix elements which are much smaller than the diagonal ones for any 4>h and a 
values (as those enter via arguments of a cosine function). We have found that for any R^^ > 400 GeV, 
if e ~ l/2i?, then all possible combinations of (ph and a yield realistic spectra. Naturally, these possible 
combinations of (j)h and a include a = as well. For possible lower values of e, which may not correspond 
to TO ^ e, a realistic spectra requires the second scenario (ph ^ + 2qTT (as already discussed). Note 
that for this case also, a = is a natural possibility. Thus the important conclusion that we made from 
these observations was that a non-zero brane-shift is not essential for a realistic mass spectrum. 

Let us compare our observations with those of Ref.[ll]. The authors of reference [11] found that when 
e = 0, (/>/j = and a = 0, the mass spectra consists of massive pairwise degenerate neutrino states without 
any massless state. They also found that when e = 1/2R, 0^ = and a = 0, there is an almost massless 
state, and the remaining spectra consists of pairwise degenerate massive neutrino states. We have found 
the same observation for these two points. 

They further found that when e ^ and ^ l/2i?, (j)h = and a = 0, there is no massless state and the 
mass spectra is hierarchical. We indeed find an almost massless neutrino state in this case too, the overall 
spectrum being hierarchical. 

They also found that unless e = l/2i?, = it/A and a — 7ri?/2, the mass spectrum consists of massive 
non-degenerate KK neutrinos. We find indeed that for e = 1/2R, (ph — 7r/4 and a = ■kR/2, the high 
neutrinos states are pairwise degenerate. 

However, we found that for several other combinations of these three parameters one can also get a realistic 
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Figure 2: Plot of m^^^/e versus e (in GeV). 

mass spectrum. For example, we mentioned before that for e = 1/2R, we always get the spectra consisting 
of an almost massless state and other states pairwise degenerate, irrespective of the numerical values of 
4>fi and a. Our larger allowed domain of parameter space is due to the additional de-correlation condition 
required in [11], and not considered here. In particular, a non-zero brane-shift parameter is found to be 
necessary in [11] in order to account for this complete de-correlation of the effective Majorana mass term 
(to) and the scale of light neutrino mass (the motivation being to have (to) in a potentially observable 
range). In contrast, we find the brane-shift to be not essential for a realistic mass spectrum. 

5. The quantity to*^"^: We finish this section by determining which of the two possible ED mechanisms 
responsible for the suppression of m^-^ (discussed in Sec3) plays the dominant role (for different values of 
R~^). Our discussion is based on the plot of rrS'^^ /e versus e in Fig. 2. One first observes clearly from 
the figure that systematically to'-"-' < 10~^£ which means that the effective ED see-saw mechanism plays 
always a significant role in reducing TOj/^ (see Eq.(22)), as we have already noticed. One see typically 
that the m,j^ suppression, relative to the electroweak symmetry breaking scale {v ~ 10^ GeV), can be 
due at 50% from the effective ED see-saw mechanism and at 50% from the higher dimensional mechanism 
based on wave function overlap, a situation corresponding to e 10^ GeV and to^"^ ~ 10~^ GeV (such 
that m^^^/v ^ to'-'^^/e ^ 10~^). The suppression can even come purely from the effective ED see-saw 
mechanism, namely to'"'/^ ~ 10~^^ with to^^^ v ^ 10^ GeV. Thus the fact that m'^^^/e < 10^^ means 
that the ED see-saw mechanism is systematically the dominant one. 

This can interpreted in the following terms. All the conditions that we have imposed resulted in a lower 
bound on (see above)^. Now, as R^^ increases, the extra compact bulk space decreases, so that the 
wavefunction overlap factor increases, making the suppression of m^^ from the geometrical mechanism 
less effective. Hence, the different constraints (including the SNO one) limit the wave function overlap 
mechanism effect. Nevertheless, thanks to the ED see-saw mechanism, a realistic spectrum can still be 
generated. 

4.2 Results for the Ug > 1 scenario 

We consider now another scenario where the right handed neutrino propagates only along one extra dimension 
whereas gravity propagates in the whole bulk, the number ng and sizes of extra dimensions being adjusted (with 
Eq.(16) when all of them have the same size 2nR) such that Mp is as low as ~ C(l) TeV in order to address 

^In particular, the SNO bound constrains the mixing angles between active and sterile neutrinos, by this way forcing typically 
the KK masses (oc to increase. 



10 




Figure 3: Diagrams for the CP-asymmetry decays, at tree (a) and loop (b)-(c) level, of the second neutrino 
eigenstate V2 into the charged lepton and the charged Higgs boson. 

the gauge hierarchy problem. Recall that the solution to the gauge hierarchy problem within the ADD model 
for rig = 1 is excluded by experimental arguments relative to the Newton law modification'^. 

The radius, as found previously, is bounded from above by R"^ > 400 GeV, since R^^ > 2|e|, e w rrii.^ 
and m^^ > 200 GeV. As a consequence, m^°^ ~ v/^/MJR > 10^ GeV ^ as Mf ~ 0(1) TeV, so that there 
is no significant suppression in m^*^^ (and hence in m^-^), relative to the electroweak scale v, coming from 
the ED overlap mechanism. Thus the neutrino mass suppression must originate principally from the ED 
see-saw mechanism. For that purpose, one should have m'^^^/e ~ 10~^^ (see Eq.(22)) which together with 
777,^°) ~ v/y/MpR and |e| < l/2i? leads to R"^ > 10^^ GeV: an amount much larger than the fundamental 
scale Mf ^ 0(1) TeV. Thus, for the ng > 1 scenario, generating a realistic neutrino spectrum and addressing 
the gauge hierarchy issue simultaneously, is only possible if there exists a compactification scale 1/R extremely 
different from the fundamental scale Mp (hence another hierarchy has to be introduced). Besides, such a scenario 
is automatically in a non-predictable regime, since the condition Afmax < Mp translates into 1/R TeV. 



5 Leptogenesis 

One interesting question to address is the possibility of having successful leptogenesis in these ED scenarios. 
As has been argued in [6] , the infinite KK tower of neutrinos associated with the 5-dimensional sterile neutrino 
may act as an infinite series of CP-violating resonators. Thus, we probed the ED scenarios to see whether they 
can generate sufficient lepton asymmetry. 

5.1 Additional Constraints and Expectations 

In the previous section we saw that the lightest neutrino vi in the generated spectrum fulfills the absolute scale 
bound. We also constrained the next to lightest neutrino 1^2 to be heavier than the Higgs mass, so that it may 
decay at a temperature corresponding to its mass and produce lepton asymmetry [13]. The out-of equilibrium 

^In our case for the Ug = I scenario, the cutoff scale was set around 10^"^ GeV and more. This model, for sure, did not address 
the gauge hierarchy issue. 

''Unless of course iph ~ ''^/'^ + 2g7r, but in order to decrease significantly the order of magnitude of this lower bound on m^o), <t>h 
should be extremely close to n/2 + 2qn which has a priori no theoretical justification. 
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condition may be written as the following condition on the decay width of f2, 



Sir 



<H{T = m.,), K.=[X'^\]. 



(26) 



where H is the Hubble expansion rate at T 



1 ^'^d is modified Yukawa coupling for the neutrino 

,0) 



mass eigenstate uj. The initial Yukawa couplings, A', can be read from Eq.(17): m^-'-' = XjV/y^. When the 
complex symmetric mass matrix in Eq.(16) is diagonalized, the weak states 'I'l, of Eq.(14) transform to the 
mass eigenstates vj/^ass fQ^Q-^ying -qi^ = y^j^mass^ where V is the diagonalizing matrix. Thus the initial Yukawa 
couplings for right handed neutrinos A^ also transform into the modified Yukawa couplings Xj via the V matrix. 
We note the Yukawa couplings Xi instead of A^i, where the index 1 is fixed as it corresponds to the unique 
left handed SM lepton (we remind that we work under the simplification assumption of one family) which 
moreover does not develop a KK tower. The condition in Eq.(26) is the first filter to ensure a successful thermal 
leptogenesis. The other necessary condition is to generate enough CP asymmetry, the asymmetry being defined 

by, 



r(i^2 ^ e^) - r(z/2 ^ l4>) 



r(i^2 



r(j/2 



(27) 



This asymmetry gets a non-zero value due to the interference between the tree level diagram of Fig. 3a and 
one-loop diagrams of Fig. 3b and 3c. The asymmetry is given by. 



j¥2 



(28) 



where / is the loop factor [24], given by 



f{x) = s/x 



l_(l + a;)lnl±^ + 



1 -X 



(29) 



in the case of a hierarchical neutrino spectrum. Here i>2, the next-to-lightest neutrino, is the external (i.e. 
decaying) mode and Vj the exchanged eigenstate (see Fig. 3). The amount of CP asymmetry generated can be 
magnified by the loop function and/or the Yukawa couplings. However, due to the large hierarchy between the 
lightest first state i>i and the next-to- lightest state 1^2, much enhancement from the Yukawa couplings cannot 
be expected, as we will discuss more precisely in the next subsection. The other enhancement may then come 
from the loop function via resonances. The discussion on the characteristics of the generated KK neutrino 
spectrum in the previous section clearly shows that we are indeed often in this favorable resonance situation 
{'^ e ^ since the neutrino eigenstates are almost pairwise degenerate over a wide range of the obtained 

parameter space. In this situation, the self energy contribution to the loop factor, which is the last term in 
Eq.(29), is modified to [25] 



(ml. 



TO^ )m. 



1^2 Vj 



(30) 



where F^^ ~ Kj^rrivJ^i^ is the decay width of state Uj. 



5.2 Results and Discussion 

As before, we have scanned the total parameter space with the additional constraints mentioned above concern- 
ing leptogenesis, concentrating on the scenario with Ug = 1. We checked that within the previously obtained 
parameter space, the first test to ensure leptogenesis (represented by Eq.(26)), is satisfied. 
As shows Eq.(28), the source of a non-zero CP asymmetry is the difference of phases of the Yukawa couplings. 
Since we are working with one generation of neutrinos, there is a single (common) phase associated with the 
Yukawa couplings A' (the one associated with the SM neutrino). However, the final (modified) Yukawa couplings 
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Xj may still be different for different neutrino mass eigenstates , as V itself is a complex matrix. The resulting 
difference of phases, generated in the process of diagonalization of the mass matrix as described above, leads to 
a non-zero value of Im[A^A]2j- 

However, although non-zero, Irn[A''^A]|^- is small due to the large hierarchies among the elements of the mass 
matrix, which in turn is essential to generate a realistic neutrino spectrum with small active-sterile mixing 
angles. We find numerically that the resonance enhancement in e^p from the loop function cannot compensate 
this suppression from the Yukawa couplings in order to ultimately generate a sufficient CP asymmetry for a 
successful leptogenesis: e^p > 10~^ [26] (assuming the usual orders of magnitude for the washout factors). Thus, 
although the pairwisc quasi-degenerate KK neutrino spectrum for one-generation of SM neutrino is favorable 
for yielding a good value of e^p , but the suppression of the same from the Yukawa couplings turns out to be 
too dominant. 

A possible alternative for a potentially successful leptogenesis could be to probe a more realistic model with 
three generations of SM neutrinos on the 3-brane, with three extra right handed neutrinos in the bulk. The mass 
matrix in Eq.(16) in that case would have a similar structure with hi and h2 replaced by 3 x 3 matrices in the 
flavor space. The masses of the n}^ neutrino states would remain approximately the same for each generation 
and the differences in the masses of the three light neutrinos would only come from the Yukawa couplings. This 
initial difference in the Yukawa couplings could be enhanced further during the process of diagonalization of 
the mass matrix. This is in contrast to the one-generation case where we had to start from a single Yukawa 
coupling, and the process of matrix diagonalization was the only source to create a difference in the final Yukawa 
coupling phases. Therefore, we could expect a much larger contribution to the CP asymmetry from the Im[A^A] 
factor, compared to the one-generation case. 



6 Conclusions 

In this work, we have considered an L-breaking extra dimensional scenario, with the SM particles localized on a 3- 
brane and an additional right handed neutrino propagating along an extra dimension, providing new mechanisms 
for suppressing the neutrino mass. We probed the possibilities of this model to account simultaneously for a 
realistic neutrino mass spectrum and a sufficient lepton asymmetry (from the decays of the KK neutrinos). We 
considered two variants of the above scenario. 

The first is when this extra dimension is the only available one {rig = 1). Then there exist fundamental 
parameters which give rise to a neutrino spectrum respecting the experimental constraints, on neutrino masses 
and active-sterile neutrino mixings (like the SNO bounds). We have found that in this framework, the neutrino 
mass suppression originates dominantly from the ED see-saw mechanism and partially from the ED wave 
function overlap effect. In fact, the SNO data limit the effect of the wave function overlap but a sufficiently 
small neutrino mass can be generated thanks to the ED see-saw mechanism. For certain parameter domains, 
the spectra in this scenario consists of pairwise quasi-degenerate heavy (KK) neutrinos, together with a light 
neutrino. We observed that under the constraints on neutrino masses that we have imposed, a "brane shifted" 
framework of the extra dimensional model is not essential for generating a realistic mass spectra. Also, the 
inverse radius of compactification has to be pushed to ~ 400 GeV (or more). This in turn results in very tiny 
phase differences among the modified Yukawa couplings (after the process of mass matrix diagonalization) so 
that it is not possible to create a CP asymmetry large enough to insure a successful leptogenesis. Indeed, such 
a tininess is not sufficiently compensated by the enhancement from the resonant loop function (related to the 
decay of the KK neutrinos) which is due to the pairwise quasi-degeneracy property of the spectrum. 

In the se(;on(l variant where gravity can propagate in more extra dimensions (rig > 1), the fundamental gravity 
scale is reduced down to the TeV scale so that the gauge hierarchy problem is addressed through the ADD 
approach. We have found that in this framework, the neutrino mass suppression originates also primarily from 
the ED see-saw mechanism and the compactification scale is of an order of magnitude much larger than the 
fundamental gravity scale. 
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A: Approximate analytical estimation of KK-neutrino masses 



Here we intend to show analytically that when the eigenvalues (5„ with n > 2) of mass matrix (16) are taken 
equal to the diagonal elements, their correction are much smaller than the eigenvalues themselves. For that 
purpose, let us start with the 4x4 form of Eq.(16), 
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(31) 



All the obtained sets of parameters (see discussion on Fig.l) lead to a significant hierarchy in m/e (case (1) of 
Eq.(23)). Hence the m^") values are systematically much smaller than e, from Eq.(17), namely: m^"-* <C £. We 
can thus rewrite matrix in Eq.(31) as = M.D+ Mi, where Md '^s, the diagonal matrix constructed from M., 
and M.I is a perturbation matrix with no diagonal entries given as. 
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(32) 



The zeroth order eigenvalues £i of M can be written as, f i = 0, ^2 = £, ^3 = £ — ^4 = £ + l/R- Since Mi 
has all diagonal elements equal to zero, the lowest order correction to the eigenvalues reads as. 



S£n= Yl 



{<i^n\Ml\i,n- >f 

£ — £ I 



(33) 



where the V's are the zeroth order eigenvectors (considering the non-degenerate case). Thus, for our case we 
have from the above relations, 



6£i = 



(mW)^ ^ (m(^))^ ^ (m(-i))2 



(34) 



l/R -e+l/R 

(35) 

5£„ - (™("')Vfn, (36) 
where m*^") is given by Eq.(17) and n > 1. Since m <C £, the first two of the above relations lead to 

2£ r 

£ 



{llRY-{ey 



S£i < 
5£2 < rr?le. 

These inequalities can immediately be generalized to the N x N mass-matrix as 



(37) 
(38) 



5£i < 



\N-1)/2 
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(n/i?)2 - (£)2 



5£2 < m^/e. 



(39) 

(40) 
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For the maximum value of e (= 1/2R), the first term on the right hand side of Eq.(39) is ~ 2R for large N, so 
that 5£\ ^ 0. For smaller values of e, SSi would be dominated by the second term on right hand side. For all 
the other eigenvalues with n > 2, one can derive £n > Thus, in general the numerical value of S£n for any n 
verifies 

\5£„\ < m^/e. (41) 

So let us study the m'^/e value. For < 10^ GeV, e > 200 GeV (sec Fig.l) and m < IQ-^ GeV (c./. Eq.(18) 
and Eq.(5) in case of the scenario with ng = 1). Hence, one has rn^/e < 10~^ GeV. Now for > 10^ 
GeV, e ~ R-^/2 and m ~ 10"'' {R-^/lGeVy/^ GeV, so that m'^/e ~ IQ-^ (i?-VlGeV)-i/3 GeV and thus 
rn^/s ^ 10~^^ GeV. Therefore, from Eq.(41), we conclude that one has systematically \6£n\ << £n for each 
n > 2, since £n>2 ^ £ ^ 200 GeV. A similar analysis could be performed for the case of the degenerate mass 
spectra (then we would have to start from the corresponding relation of Eq.(33) for degenerate eigenvalues). 



15 



References 



[1] N. Arkani-Haincd, S. Dimopoulos and G. Dvali, Phys. Lett. B 429 (1998) 263; I. Antoniadis, N. Arkani- 
Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B 436 (1998) 257; K.R. Dienes, E. Dudas and 
T. Gherghetta, Phys. Lett. B 436 (1998) 55; Nucl. Phys. B 537 (1999) 47. 

[2] L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 3370. 

[3] L Antoniadis, Phys. Lett. B 246 (1990) 377; J.D. Lykken, Phys. Rev. D 54 (1996) 3693; P. Hofava and 
E. Witten, Nucl. Phys. B 460 (1996) 506; Nucl. Phys. B 475 (1996) 94. 

[4] K.R. Dienes, E. Dudas and T. Gherghetta, Nucl. Phys. B557 (1999) 25. 

[5] N. Arkani-Hamcd, S. Dimopoulos, G. DvaU and J. March-Russell, Phys. Rev. D 65 (2002) 024032. 

[6] A. Pilaftsis, Phys. Rev. D 60 (1999) 105023. 

[7] A. loannisian and A. Pilaftsis, Phys. Rev. D62 (2000) 066001. 

[8] M. Maltoni et al, New. J. Phys. 6 (2004) 122. 

[9] For example, see: G.R. DvaU and A.Y. Smirnov, Nucl. Phys. B 563 (1999) 63; R.N. Mohapatra, S. Nandi 
and A.Perez-Lorenzana, Phys. Lett. B 466 (1999) 115; A. Lukas, P. Ramond, A. Romanino and G.G. Ross, 
Phys. Lett. B 495 (2000) 136; R. Barbieri, P. CreminelU and A. Strumia, Nucl. Phys. B 585 (2000) 28; 
R.N. Mohapatra and A.Perez-Lorenzana, Nucl. Phys. B 593 (2001) 451; A. Lukas, P. Ramond, A. Romanino 
and G.G. Ross, JHEP 0104 (2001) 010; K.R. Dienes and I. Sarcevic, Phys. Lett. B 500 (2001) 133; 
N. Cosme, J.-M. Frere, Y. Gouverneur, F.S. Ling, D. Monderen and V.Van Elewyck, Phys. Rev. D 63 (2001) 
113018; D.O. Caldwell, R.N. Mohapatra and S.J. Yellin, Phys. Rev. Lett. 87 (2001) 041601; J. Maalampi, 
V. Sipilainen and I. Vilja, Phys. Lett. B 512 (2001) 91; C. S. Lam and J. N. Ng, Phys. Rev. D 64 (2001) 
113006; A.S. Dighc and A.S. Joshipura, Phys. Rev. D 64 (2001) 073012; A. Nicolaidis and D.T. Papadamou, 
Phys. Rev. D 66 (2002) 013005; J.-M. Frere, M.V. Libanov and S.V. Troitsky, JHEP 0111 (2001) 025; 
J.L. Hewett, P. Roy and S. Roy, Phys. Rev. D 70 (2004) 051903; N. Arkani-Hamed and S. Dimopoulos, 
Phys. Rev. D 65 (2002) 052003; T. Gherghetta, Phys. Rev. Lett. 92 (2004) 161601; N. Haba, JHEP 0605 
(2006) 030; M.-T. Eisele and N. Haba, Phys. Rev. D 74 (2006) 073007. 

[10] A. Lukas, P. Ramond, A. Romanino and G. G. Ross, JHEP 0104 (2001) 010. 

[11] G. Bhattacharyya, H.V. Klapdor-Kleingrothaus, H. Pas and A. Pilaftsis, Phys. Rev. D 67 (2003) 113001. 

[12] H. Davoudiasl, P. Langacker and M. Perelstein, Phys. Rev. D 65 (2002) 105015. 

[13] M. Fukugita and T. Yanagida, Phys. Lett. B 174 (1986) 45. 

[14] For example, see, E.G. Gimon and J. Polchinski, Phys. Rev. D 54 (1996) 1667. 

[15] Y. Grossman and M. Neubert, Phys. Lett. B 474 (2000) 361; G. Moreau, Eur. Phys. J. C 40 (2005) 539. 

[16] M. CireUi, G. Marandella, A. Strumia and F. Vissani, Nucl. Phys. B 708 (2005) 215. 

[17] P. Aliani, V. Antonelli, M. Picariello and E. Torrente-Lujan, arXiv:hcp-ph/0309156. 

[18] A. Bandyopadhyay, S. Choubey, S. Goswami, S.T. Petcov and D.P. Roy, Phys. Lett. B 583 (2004) 134. 

[19] P. Minkowski, Phys. Lett. B 67 (1977) 421; M. GcU-Mann, P. Ramond and R. Slansky, Proceedings of 
the Supergravity Stony Brook Workshop^ New York 1979, eds. P. Van Nieuwenhuizen and D. Frcedman; 
T. Yanagida, Proceedings of the Workshop on Unified Theories and Baryon Number in the Universe, 
Tsukuba, Japan 1979, ed.s A. Sawada and A. Sugamoto; R.N. Mohapatra, G. Senjanovic, Phys. Rev. 
Lett. 44 (1980) 912. 

[20] SNO Collaboration, Phys. Rev. Lett. 87 (2001) 071301. 



16 



[21] Supcr-Kainiokande Collaboration, Phys. Rev. Lett. 86 (2001) 5651. 

[22] T. K. Kuo and J. Pantaleone, Rev. Mod. Phys. 61 (1989) 937. 

[23] See, e.g. N.F. Bell, FERMILAB-FN-0785-A, T.M. Ito and G. Prezeau, Phys. Rev. Lett. 94 (2005) 161802. 

[24] L. Covi, E. Roulet and F. Vissani, Phys. Lett. B 384 (1996) 16. 

[25] A. Pilaftsis and T.E.J. Underwood, Nucl. Phys. B 692 (2004) 303. 

[26] See, e.g. G.F. Giudice, A. Notari, M. Raidal, A. Riotto and A. Strumia, Nucl. Phys. B 685, 89 (2004); 
W. Buchmuller, P. Di Bari and M. Plumacher, Annals Phys. 315 (2005) 305. 

[27] V.A. Rubakov and M.E. Shaposhnikov, Phys. Lett. B 125 (1983) 136. 

[28] N. Arkani-Hamed and M. Schmaltz, Phys. Rev. D 61 (2000) 033005; E.A. Mirabelli and M. Schmaltz, 
Phys. Rev. D 61 (2000) 113011; G. Dvah and M. Shifman, Phys. Lett. B 475 (2000) 295; G.C. Branco, 
A. de Gouvea and M.N. Rebelo, Phys. Lett. B 506 (2001) 115; M. Raidal and A. Strumia, Phys. Lett. B 
553 (2003) 72. 

[29] H.V. Klapdor-Kleingrothaus and U. Sarkar, Phys. Lett. B 541 (2002) 332. 



17 



